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Abstract 

Unsteady transonic flow calculations are 
presented for aerodynamical ly Interfering 
lifting surface configurations. Calculations 
are performed by extending the XTRAN3S (Version 
1.5) unsteady transonic small -disturbance code 
to allow the treatment of an additional lifting 
surface. The research was conducted as a 
first-step toward developing the capability to 
treat a complete flight vehicle. Grid genera- 
tion procedures for swept tapered Interfering 
lifting surface applications of XTRAN3S are 
described. Transonic calculations are presented 
for wing-tail and canard-wing configurations for 
several values of mean angle of attack. The 
effects of aerodynamic interference on transonic 
steady pressure distributions and steady and 
oscillatory spanwise lift distributions are 
demonstrated. Results due to wing, tall, or 
canard pitching motions are presented and 
discussed in detail. 


Nomenclature 

c airfoil chord 

c 4 sectional lift coefficient 

c r wing reference chord 

f function defining instantaneous position 

of lifting surface 
k reduced frequency, wc r /2U 

M freestream Mach number 

S horizontal separation distance measured 

in wing root chordlength between the two 
planform root sections 

t time nondimensionalized by freestream 

velocity and wing reference chord 
U freestream velocity 

x, y,z nondimensional cartesian coordinates In 

streamwlse, spanwise, and vertical 
directions, respectively 

y, coordinate of jth spanwise grldpolnt 

a angle of attack 

y ratio of specific heats 

aCp lifting pressure coefficient 

A v sweep angle 

"n fractional semi span 

5,n,c transformed nondimensional coordinates 

in x, y, and z directions, respectively 
<p disturbance velocity potential 

<*> angular frequency 

Subscripts 

C canard 

LE leading edge 

T tail 

TE trailing edge 

W wing 

0 mean value 

1 dynamic value 


Introduction 

Computational methods employing linear 
theory for predicting unsteady flowflelds about 
aerodynamical ly interfering lifting surfaces and 
multiple-component configurations have been 
developed over the past fifteen years (Refs. 
1-4, for example). These methods are extensions 
of linearized theory for single lifting surfaces 
to more complicated configurations such as two 
lifting surfaces In tandem, various horizontal 
stabilizer/vertical fin arrangements, and wing- 
body combinations. More recently, computer 
programs such as that reported in Ref. 5 have 
been developed for the analysis of complete 
flight vehicles. The underlying linear theory 
on which these methods are based, though, 
restricts the applications to either subsonic or 
supersonic flows. In the transonic regime, 
computational methods for predicting flowflelds 
about oscillating multiple lifting surface 
configurations are currently being developed. 

Steady transonic flowflelds about isolated 
wings have been calculated by Boppe b using a 
finite-difference relaxation scheme with a grid 
embedding procedure. The isolated wing formula- 
tion was extended to permit the analysis of two 
interfering wing surfaces. Steady pressure 
distributions were presented to illustrate 
canard-wing interference at transonic speeds. 
Shankar and Malmuth 7 presented computations for 
a few canard-wing configurations using transonic 
small-disturbance theory. In Ref. 7, results 
were obtained by placing the two surfaces in a 
single sheared fine grid system that is embedded 
in a global cartesian crude grid. Shankar and 
Goebel b developed a local numerical mapping 
procedure where the leading and trailing edges 
of the two surfaces are treated as constant 
coordinate lines in the computational domain. 
The computations have also recently been 
extended to include dihedral and fuselage 
effects. 9 

Unsteady transonic results for two-dimen- 
sional canard-wing configurations were reported 
by the author in Ref. 10. The aerodynamic cal- 
culations were performed by extending the 
XTRAN2L 11 transonic small -disturbance code to 
include an additional airfoil. Results for a 
range of canard-wing separation distances 
demonstrated the effects of aerodynamic Inter- 
ference on transonic unsteady airloads, aero- 
elastic stability, and flutter of oscillating 
airfoils. 

This paper presents three-dimensional 
transonic unsteady aerodynamic results from the 
XTRAN3S 12 transonic small -disturbance (TSD) code 
for oscillating interfering lifting surface 
configurations. The objectives of the study 
were to: (1) develop and demonstrate an 
unsteady transonic computational capability for 
Interfering lifting surface geometries; and (2) 
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investigate the effects of aerodynamic 
interference on transonic pressures and forces. 
In this study, the calculations are performed by 
extending the XTRAN3S code to allow the treat- 
ment of an additional lifting surface. The 
interfering lifting surface computational 
capability now permits the assessment of aero- 
dynamic Interference effects on transonic 
unsteady airloads and flutter characteristics of 
closely-coupled wing-tall and canard-wing 
systems. 


Computational Procedures 

In this section, the XTRAN3S code and 
coordinate transformation details are 
described. Grid generation procedures for 
Interfering lifting surface applications are 
also given. 

XTRAN3S Code 

Program Description . - The XTRAN3S code 
provides a time-accurate solution to the 
nonlinear, small -disturbance, potential equation 
for three-dimensional transonic flow. The code 
can be used to calculate steady and unsteady 
transonic flowflelds about isolated planar wings 
including aeroelastlc deformation effects. The 
program Is capable of treating either forced 
.harmonic or aeroelastlc transient type motions. 
A time-accurate alternating-direction implicit 
(ADI) finite-difference scheme Is used to solve 
the modified transonic small -disturbance 
equation 

M 2 U t + 2* x ) t = [(1 - M 2 )$ x + A* 2 + B4 > 2 ] x 

( 1 ) 

+ ( <f>y + C'VyJy + 

Several choices are available for the 
coefficients A, B, and C, depending upon the 
assumptions used In deriving the TSD equation. 12 
Briefly, the coefficients are referred to as 


“NASA Ames" coefficients when defined as 


A » - \ (y + 1)M 2 

(2a) 

B = \ (y - 3)M 2 

(2b) 

C = -(y - 1)M 2 

(2c) 

and are referred to as "NLR" coefficients 
defined as 

when 

A = - \ [3-(2-y)M 2 JM 2 

(3a) 

B = - | M 2 

(3b) 

C = - M 2 

(3c) 

The "classical" coefficients are given by 


A = - \ (y+1)M 2 

(4a) 

B = 0 

(4b) 

O 

II 

O 

(4c) 


and 'finally the coefficients for the linear 
equation are 

A = B = C - 0 (5) 

The subsonic flat plate results presented herein 
are computed using the linear equation 
coefficients; the transonic results are computed 
using the NASA Ames coefficients. For wing-tall 
and canard-wing applications, the ADI solution 
procedure of XTRAN3S has been extended to allow 
the treatment of an additional lifting surface. 
The modifications have been made to Version 
1.5 13 of the XTRAN3S code. The program now is 
capable of computing unsteady transonic 
flowflelds about relatively general Interfering 
lifting surface configurations. The present 
program coding, however, does not allow for 
overlapping or co-planar configurations. 

Boundary conditions Imposed upon the 
flowfleld are 


Far upstream: 

♦ = 0 

(6a) 

Far downstream: 

*x + \ = 0 

(6b) 

Far above and below: 

= 0 

(6c) 

Far spanwise: 

*y "0 

(6d) 

Root reflection pl^ne 

♦y = 0 

(6e) 

Trailing wakes: 

[♦ z ] - o 

(of) 


[$ x + * t ] = o 

(6g) 

where [ ] indicates 

the jump In the 

Indicated 


quantity across the wake. The flow tangency 
boundary condition is 


which is applied on the respective mean planes 
of the lifting surfaces. 

Coordinate Transformation. - The 
finite-difference grids In both the physical and 
computational domains are contained within 
rectangular regions and conform to the two 
lifting surfaces. Regions In the physical 
domain, such as the swept and/or tapered lifting 
surface planforms and their outboard extensions, 
are mapped Into rectangular regions In the 
computational domain using the shearing 
transformation 


e = 5(x.y) 

(8a) 

n = y 

(8b) 

5 = z 

(8c) 


The coordinate transformation (Eqs. (8)) is 
based on the single surface transformation 
equations described in detail by Bennett, et 
al. 14 For Interfering lifting surface 
configurations, the transformation Is applied 
independently In five separate regions of the 
grid as Illustrated In Fig. 1. The five regions 
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Fig. 1 Regions in physical domain which 

transform to rectangular regions In 
computational domain. 


are defined as: (1) upstream of the "forward" 
planform; (2) along the forward planform and 
outboard extension; (3) between the two 
planforms; (4) along the "aft" planform and 
outboard extension; and (5) downstream of the 
aft planform. 


Grid Generation Procedures 

Sectional Grid Treatment. - The sectional 
finite-difference grid near a simple wing-tail 
configuration is shown in Fig. 2. Planar wakes 
are assumed for both lifting surfaces. The grid 
is based on the XTRAN3S isolated wing default 
grid described in Ref. 15 and has 93 and 46 
points in the horizontal and vertical 
directions, respectively, for the configuration 
shown. Thirty equidistantly spaced points are 
distributed along both the wing and tall 
chords. One additional point is added near the 
leading edge of each surface for better 
resolution. Grid stretching to the farfield In 
both the x and z-di recti ons is identical to that 
described in Ref. 15. In the streamwise 
direction, for example, gridpoints are stretched 
twenty chordlengths to the farfield with eleven 
and ten points extending to upstream and 
downstream boundaries, respectively. Above the 
tail and below the wing, twenty gridpoints are 
stretched twenty-five chordlengths to the upper 
and lower farfield boundaries. Between the two 
surfaces, vertical grid lines are distributed 
equidistantly between the first gridpoint 
downstream of the wing trailing edge and the 
tail leading edge. Horizontal grid lines are 
distributed equidistantly between the wing and 
tail centerlines with one additional line placed 
symmetrically below the tail centerline and 
another additional line placed symmetrically 
above the wing centerline. Horizontal and 
vertical separation distances are 50% and 40% 
chord, respectively, for the configuration shown 
here. For different separation distances or 
other configurations, points may be added or 
removed between the two lifting surfaces. 

Planform Grid Treatment. - Figure 3 shows 
the planform finite-difference grid near a 
simple wing-tall configuration consisting of two 
Identical swept planforms. The grid is based on 
the single surface grid mapping procedures 
described in Ref. 14 which were extended in this 


study to include an additional planform. A 
cosine distribution of spanwise gridpoints is 
used in this example along the lifting surfaces 
and in the outboard region to concentrate points 
near the tips. Outboard of the lifting 
surfaces, the planforms are extended to the 
farfield boundary using the outboard grid 
mapping procedure of Bennett, et al. 11 * In this 
procedure, the leading and trailing edges are 
smoothly extended using cubic polynomials that 
match leading and trailing edge slopes and 
Intersect the far spanwise boundary 
perpendicularly. The spanwise grid lines are 
mapped into the computational domain where they 
are treated as constant coordinate lines, 
similar to the local numerical mapping procedure 
of Ref. 8. Details of the grid mapping 
procedure and the equations defining the 
gridpoints are given in Ref. 14 for single 
lifting surface applications of XTRAN3S. 
Planform grids for the interfering lifting 
surface configurations presented below 
graphically illustrate the outboard grid mapping 
procedure. 

In Initial applications to both single and 
interfering lifting surface cases, numerical 
instability problems at the root were 
encountered. A representative example for a 
forward-swept wing configuration is shown in 
Fig. 4(a). Oscillations in the solution develop 
at the root as shown for t = 1.8 (120 time 
steps). These oscillations grow large with time 
and subsequently lead to program failure at t = 
2.43 (162 time steps) for this example. A study 
was undertaken to investigate this instability 



Fig. 2 XTRAN3S sectional finite-difference grid 
near wing-tail configuration of Ref. 4. 



Fig. 3 XTRAN3S planform finite-difference grid 
near wing-tail configuration of Ref. 4. 


3 



and to seek solutions to eliminate the problem. 
The source of the Instability was traced to the 
original treatment of the symmetry condition and 
grid near the root reflection plane for certain 
configurations and flow conditions. The 
symmetry condition (Eq. 6(e)) Is expressed In 
the computational domain as 

* 5 Sy + * n - 0 (9) 

where Cy Is the spanwlse grid metric along the 
reflection plane. In XTRAN3S, a “dumny" plane 
adjacent to the reflection plane provides 
grldpoints for a central -difference 
approximation of <p n in Eq. (9). The locations 
of these grldpoints in the original code are 
determined by extrapolating the grid lines In 
the physical domain across the reflection plane 
to the dummy plane as shown in Fig. 4(a). A new 
treatment of the grid was Investigated such that 
these gridpolnts are positioned symmetrically 
about the reflection plane as shown In the left 
part of Fig. 4(b). In this new grid treatment, 
the metric Cy is discontinuous across the 
reflection plane and Is therefore not defined. 
Calculations are performed by setting equal 
to zero, which Is equivalent to taking an 
average of the values of 5y on either side of 



(a) original grid treatment at wing 
root for forward-swept wing and 
steady pressure distribution at 
t = 1.8. 



(b) synmetrlc grid treatment at wing 
root for forward-swept wing and 
steady pressure distribution at 
t = 24.0. 


the reflection plane, and the symmetry condition 
becomes simply 

. - 0 (10) 

The new symmetric treatment of the grid 
eliminates the Instability at the root and 
converged steady pressure results were obtained 
as shown In the right part of Fig. 4(b). For 
this forward-swept wing configuration, the 
program now predicts a steady pressure 
distribution with a relatively strong shock wave 
at the root. Moving outboard the shock weakens 
In strength and disappears toward the tip. 
Further Investigation Is required to assess the 
accuracy of these results. 

The symmetric grid treatment can also be 
successfully applied to aft-swept configurations 
such as that presented in Fig. 4(c). The 
resulting steady presure distribution Is smooth 
at the root and Indicates the presence of a mild 
shock that Is more pronounced at the tip. For 
highly-swept configurations, though, a non- 
physical pressure expansion may occur near the 
wing apex due to the sharp corner Introduced by 
the symmetric grid. This nonphysical pressure 
expansion occurred for the configuration of 



(c) symmetric grid treatment at wing 
root for aft-swept wing and steady 
pressure distribution at t = 24.0. 



(d) symmetric grid treatment with 
parabolic smoothing at wing root 
and steady pressure distribution at 
t * 24.0. 


resulting steady 


Fig. 4 Treatment of planform grid across reflection plane and 
pressure distributions at M * 0.9 and ct„ * 0°; 
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Fig. 4(c) when a cosine distribution of 

spanwise gridpoints was used. An Inboard grid 
mapping procedure was subsequently developed to 
eliminate the corner by smoothing the grid 
across the reflection plane. In this procedure, 
the leading and trailing edges of the planforms 
near the root sections are replaced by 

second-order polynomials of the form 

x LE (y) = \ E y 2 + B LE y + c LE (Ua) 

x TE (y) = A TE y 2 + B TE y + C TE (lib) 

The coefficients in Eqs. (11) are determined by 
requiring that the polynomials Intersect the 
reflection plane perpendicularly and match the 
locations and slopes of the leading and trailing 
edges at the first spanwise gridpoints outboard 
of the root section (y = y 3 ). The inboard grid 
mapping procedure with parabolic smoothing 
therefore replaces the leading and trailing 
edges at the root (y = y 2 = 0) by 

x LE (°) = x LE (y 3 ) - -|y 3 tan_1 (7 - a le) ( 12a ) 

x TE (°) = x TE (y 3 ) - |y 3 tan' 1 (| - a te ) (12b) 

which effectively unsweeps the grid near the 
reflection plane as shown in Fig. 4(d). This 
grid variation produces metrics Cy that are 
continuous and, along, the reflection plane, are 
uniquely defined to be zero. The resulting 
steady pressure distribution using this grid 
treatment is also shown in Fig. 4(d). These 
steady pressures are smooth at the root and are 
similar to those shown in Fig. 4(c) (calculated 
using no smoothing along the root). A 
comparison of steady pressure distributions 
computed using three different grid variations 
is shown in Fig. 5 for the root section of the 
aft-swept wing. The symmetric treatment of the 
grid with parabolic smoothing produces a steady 
pressure distribution that is more forward 
loaded in comparison with the original and 
symmetric (without root smoothing) grids. This 
effect is attenuated In the spanwise direction 
and is negligible at the tip. Further work is 
needed to determine the accuracy of results 



Fig. 5 Comparison of steady pressure 

distributions for aft-swept wing root 
section computed using three grid 
variations across reflection plane. 


computed using the symmetric grid with root 
smoothing by comparison with experiment. 

Planform finite-difference grids for the 
wing-tail and canard-wing configurations 
considered herein, are shown in Figs. 6(a) and 
6(b), respectively. These grids clearly 
illustrate the Inboard root smoothing and the 
outboard grid mapping procedures for interfering 
lifting surface configurations. The outboard 
grid mapping of Ref. 14 is applicable for most 
wing-tail and canard-wing configurations. For 
geometries consisting of a forward-swept wing 
with a closely-coupled aft-swept canard, 
however, the outboard extensions of the two 
planforms may overlap. The mapping has 
subsequently been modified to resolve this 
problem by simply altering the locations where 
the outboard extensions Intersect the far 
spanwise boundary. The modification separates 
the two outboard extensions such that the 
overlap no longer occurs and reasonable planform 
grids for closely-coupled configurations may be 
obtained. 


Results and Discussion 

Configurations 

Results are presented for the following 
three configurations: (1) a simple wing-tail 
consisting of two identical swept planforms; (2) 
the laschka-Schmidt wing-tail 3 ; and (3) a 



(a) Laschka-Schmidt wing-tail . 



(b) canard-wing configuration. 


Fig. 6 
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XTRAN3S planform finite-difference grids 
near interfering lifting surface 
configurations (S = 0.50); 


canard-wing geometry. Subsonic unsteady 
calculations were performed for the first 
configuration for comparison with the linear 
theory results of Ref. 4. Transonic steady and 
unsteady results were obtained for the second 
and third configurations to demonstrate the 
XTRAN3S interfering lifting surface capability 
and to investigate the effects of aerodynamic 
Interference on pressure distributions and 
spanwise loading. 

Subsonic unsteady results are presented for 
the simple flat-plate wing-tall configuration 
Illustrated in Figs. 2 and 3, to assess and 
verify the modifications to the XTRAN3S code. 
This configuration was selected for comparison 
with the linear theory results reported by 
Albano, Perklnson, and Rodden. 1 * In Ref. 4, 
subsonic results were obtained for this 
configuration using a kernel function method. 
As Illustrated in Fig. 3, each planform of the 
wing-tail has a full -span aspect ratio of 2.28, 
a leading edge sweep angle of 45°, and a taper 
ratio of 1.0. 

Transonic steady and unsteady results are 
presented for the Laschka-Schmidt wing-tall and 
a canard-wing configuration. Each planform of 
these configurations has a full -span aspect 
ratio of 2.67, a leading edge sweep angle of 
45°, and a taper ratio of 0.5. The tall of the 
Laschka-Schmidt configuration is two-thirds the 
size of the wing as illustrated by the planform 
finite-difference grid of Fig. 6(a). The 
canard-wing configuration is created by simply 
placing the tail in front of the wing as 

illustrated by the planform finite-difference 
grid of Fig. 6(b). The grids contain 82 x 

20 x 44 points in the x, y, and z directions, 

respectively, for a total of 72,160 points. The 
vertical separation distance for the tail and 

canard was selected as 20% wing root chord above 
the wing. Two horizontal separation distances 
(between the planform root sections) of 25% 
(S = 0.25) and 50% (S a 0.50) wing root chord 
were considered. The NACA 0010 airfoil section 
was selected for the lifting surfaces. These 
similar wing-tail and canard-wing geometries 
were selected to study transonic Interference 
effects because the lifting surfaces are all 
geometrically the same. Differences in pressure 
distributions or spanwise loading between the 
wing and the. tail or between the canard and the 
wing are, therefore, due solely to aerodynamic 
interference. 

Steady computations also were attempted for 
the Rockwell canard-wing model of Ref. 16. Each 
of the two lifting surfaces had an exposed 
aspect ratio of approximately 2.0, a leading 
edge sweep angle of approximately 40°, and a 
taper ratio slightly greater than 0.25. The 
wing also had 5° of parabolic twist washout. 
Attempts to obtain converged steady pressure 
distributions for both the canard-wing as well 
as for the isolated wing were unsuccessful due 
to a numerical Instability that first was 
observed in the wing wake. The instability was 
attributed to the large twist of the low aspect 
ratio wing since steady solutions were 
obtainable for the wing with no twist. 


Subsonic Linear Theory Comparison 

Subsonic unsteady computations for the 
swept wing-tall configuration of Figs. 2 and 3 
were performed for simple harmonic wing pitching 
about a line perpendicular to the root at the 
wing apex. The freestream Mach number was 
M = 0.7, the mean angle of attack was oq => 0° 
for both wing and tail, and the reduced 
frequency based on root semi chord was k = 
0.262. The nondlmensional time step size was 
selected as at => 0.01872 which results In 640 
time steps per cycle of motion. 

Unsteady chordwlse lifting pressure 
distributions aC p on the wing and tail due to 
wing pitching with oscillation amplitude a t = 1° 
are shown in Fig. 7 for four span stations. 
Real and imaginary components of the lifting 
pressures are shown. These results were 
obtained using the symmetric treatment of the 
grid across the reflection plane without 

smoothing. The XTRAN3S results are in good 
general agreement with the kernel function 
method results of Ref. 4, with the largest 
differences occurring in the leading edge 
regions. Calculations performed using the 
original (extrapolated) treatment of the grid, 
however, did not predict the correct singular 
behavior near the leading edge of the root 
sections. This suggests that the symmetric grid 
treatment may be more applicable to unsteady 
calculations than the original grid treatment. 
The unsteady pressures on the motionless tail 
are Induced by the oscillatory wing pitching 
motion and represent the unsteady aerodynamic 
Interference of the wing on the tall. As 
discussed by Albano, et al., 1 * the upstream 
influence of the tail on the wing Is negligible 
for this case. 

Unsteady sectional lift coefficients on the 
wing- and tail due to wing pitching are presented 
in Fig. 8. These coefficients are plotted as 
real and imaginary components of the spanwise 
lift distribution. As shown in Fig. 8, the 
XTRAN3S calculations agree very well with the 
linear theory results of Ref. 4. The unsteady 
lift coefficients on the tail Induced by the 
wing pitching are similar in magnitude to the 
self-induced lift coefficients on the wing. 



Fig. 7 Comparison between XTRAN3S and kernel 
function method unsteady chordwlse 
lifting pressure distributions on wing 
and tail due to wing pitching at M = 0.7 
and k = 0.262. 
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. Fig. 8 Comparison between XTRAN3S 'and kernel 

function method unsteady sectional lift 
coefficients on wing and tail due to 
wing pitching at M = 0.7 and k = 0.262. 


Transonic Steady Interference Effects 

Transonic results for the Laschka-Schmldt 
wing-tail of Fig. 6(a) and the canard-wing 
configuration of Fig. 6(b) were obtained at a 
freestream Mach number of M = 0.9 and three mean 
angles of attack a 0 = 0°, 1°, and 2°. For 
steady calculations, the step size was selected 
as at = 0.015 and 1600 steps were required for 
the solution to converge. 

Steady pressure distributions on the upper 
and lower surfaces of the isolated lifting 
surface at M = 0.9 and a 0 = 2° are shown In 
Figs. 9(a) and 9(b), respectively. These 
results are presented for comparison with the 
interfering lifting surface pressure 
distributions which follow. As shown in Fig. 
9(a), there is a relatively strong shock wave on 
the upper surface located near 85% chord at the 
root. In the spanwise direction, the shock 
becomes more pronounced and Is located near 30% 
chord at the tip. Although not obvious, the 
lower surface pressures shown in Fig. 9(b) are 
supercritical. Steady pressure distributions on 
the upper and lower surfaces of the wing-tall at 
M = 0.9 and a 0 = 2° are plotted in Figs. 10(a) 
and 10(b), respectively. In general, the 
pressures on the wing are very similar to the 
isolated surface pressures shown in Fig. 9. On 
the tail, however, the shock is weakened on the 
upper surface and the pressures near the lower 
surface leading edge are Increased In magnitude 
due to the aerodynamically Interfering wing. 
Also, the pressure distributions on the upper 
and lower tail surfaces are very similar which 



(a) upper surface. 



(b) lower surface. 


Fig. 9 Steady pressure distributions on 

Isolated lifting surface at M = 0.9 and 
«o = 2°; 



(a) upper surfaces. 



(b) lower surfaces. 


Fig. 10 Steady pressure distributions on 

Laschka-Schmldt wing-tail at M = 0.9 
and otg = 2° (S = 0.50); 


7 


indicates that the tail lift is small. Pressure 
distributions on the canard-wing upper and lower 
surfaces at M = 0.9 and a 0 = 2° are presented 
in Figs. 11(a) and 11(b), respectively. In 
general, the pressures on the wing resemble 
those of the Isolated surface, although the 
lower surface values are slightly Increased near 
the leading edge. On the canard, the upper 
surface shock Is slightly weakened due to the 
presence of the wing. In contrast with the tall 
pressures of Fig. 10, the canard lifting 
pressure is of sizeable magnitude since the 
upper and lower surface pressure distributions 
are dissimilar. 

To more clearly illustrate the steady 
transonic interaction between the two lifting 
surfaces, sectional lift . coefficient 
distributions are presented In Figs. 12 and 13 
for the wing-tall and canard-wing 
configurations, respectively. In each case, the 
horizontal separation distance is S » 0.5. 
Results for the Isolated lifting surface (solid 
lines) are also plotted for comparison 
purposes. For the wing-tail configuration 
(Fig. 12), the tail Induces a weak upwash on the 
wing which Increases the wing sectional lift 
coefficient C£ M for all three angles of 
attack. Conversely, the wing Induces a strong 
downwash on the tail. which significantly reduces 
the tall sectional lift coefficient c 4 _. Small 
negative values of' c^- result .for the 
an = 0° and 1° cases. Calculations for the 
closely-coupled wing-tall configuration 
(S = 0.25) resulted In negligible changes In 
C£ M and C£y (not shown). For the 




(a) upper surfaces. 



(b) lower surfaces. 



Fig. 12 Steady sectional lift coefficients 
on Laschka-Schmidt wing-tail at 
M = 0.9 (S = 0.50). 



Fig. 11 Steady pressure distributions on 

canard-wing configuration at M = 0.9 
and oq - 2° (S = 0.50); 


Fig. 13 Steady sectional lift coefficients 
on canard-wing configuration at 
M - 0.9 (S = 0.50). 
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canard-wing configuration (Fig. 13), the weak 
upstream influence of the wing slightly lowers 
the canard sectional lift coefficient c i( . at 
all three angles of attack. At a 0 = 0°, the 
aerodynamically interfering canard Induces a 
small positive spanwise lift distribution on the 
wing. As angle of attack Is increased, the 
effect of the canard on the wing Is Increased, 
especially in the inboard region. At = 2°, 
for example, the Inboard wing lift is 
significantly lowered by the canard downwash. 
The outboard wing lift is less affected by the 
presence of the canard, as expected. For the 
closely-coupled canard-wing configuration 
(S = 0.25), c iw in the inboard region of the 
wing is further reduced with decreased 
separation distance (not shown). 


Transonic Unsteady Interference Effects 

Transonic unsteady calculations for the 
wing-tail and canard-wing configurations were 

performed for wing, tail, and canard pitching 
motions with oscillation amplitude a, = 1°. The 
freestream Mach number was M = 0.9, the mean 
angle of attack was a 0 = 2°, and the reduced 
frequency was k = 0.3. The step size was chosen 
as at = 0.01496 which results in 700 time steps 
per cycle of motion. Calculations were 

performed for two different separation distances 
S = 0.25 and S = 0.50 to investigate the effects 
of separation distance on unsteady spanwise 

loadings. These calculations' are compared with 
the isolated lifting surface pitching results to 
demonstrate the effects of unsteady aerodynamic 
interference. 

Unsteady sectional lift coefficients for 

the wing-tail configuration are presented in 
Figs. 14 and 15 for wing pitching and tail 
pitching, respectively. As shown in the upper 
part of Fig. 14, the wing sectional lift 
coefficient due to wing pitching, termed 
self-induced wing lift, shows a small increase 
in the real part in comparison with the isolated 
lifting surface solution. The harmonic motion 
of the wing produces an oscillatory flow about 
the fixed tail. The resulting induced tail 

loading due to wing pitching is plotted In the 
lower part of Fig. 14. The tail sectional lift 
coefficients are of similar magnitude to the 
self-induced wing lift and are dependent upon 
the horizontal distance separating the two 
lifting surfaces. When the separation distance 
is reduced by a factor of two, the Imaginary 
part of c* T decreases and the magnitude of 
the real part increases. These changes in the 

real and imaginary parts occur such that the 
magnitude of Cj T is uneffected and the phase 
of c^-p is decreased with reduced separation 
distance S. The decreased phase of c a . Is 
physically related to the shorter time required 
for disturbances created by the wing to 
influence the tall. For tail pitching (Fig. 
15), the tail spanwise loading shows only a 
small increase in both real and imaginary parts 
between the isolated lifting surface and 

wing-tail cases. The effect of the oscillating 
tail on the fixed wing is negligible even for 
the closely-coupled case of S = 0.25. This 
result may be attributed to the relatively long 
times associated with upstream propagating 
disturbances (from tail to wing) in comparison 
with downstream propagating disturbances (from 



Fig. 14 Unsteady sectional lift coefficients 
on Laschka-Schmidt wing-tail due to 
wing pitching at M = 0.9, oq = 2°, 
and k = 0.3. 



Fig. 15 Unsteady sectional lift coefficients 
on Laschka-Schmidt wing-tail due to 
tail pitching at M = 0.9, <xq = 2°, 
and k = 0.3. 
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wing to tall). For the freestream Mach number 
of M = 0.9 considered here, the downstream 
propagation speed Is approximately nineteen 
times faster than the upstream propagation 
speed. 

Unsteady sectional lift coefficients for 
the canard-wing configuration are presented In 
Figs. 16 and 17 for canard pitching and wing 
pitching, respectively. As shown In the upper 
part of Fig. 16, the canard sectional lift 
coefficient due to canard pitching, termed 
self-induced canard lift, shows only small 
differences between the Isolated lifting surface 
and canard-wing cases. The harmonic motion of 
the canard produces an oscillatory flow about 
the wing. The resulting Induced wing loading 
due to canard pitching Is presented In the lower 
part of Fig. 16. The wing sectional lift 
coefficients are of similar magnitude to the 
self-induced canard sectional lift coefficients 
and are dependent upon the horizontal separation 
of the two lifting surfaces. With decreased 
separation distance, the magnitude of the wing 
sectional lift coefficients Is slightly 
Increased. The effect of the canard on the wing 
is largest in the inboard region and decreases 
outboard of the canard tip. For wing pitching 
(Fig. 17), the wing sectional lift coefficients 
show only very small Increases In the real and 
Imaginary components when the canard is Included 
In the transonic flowfield. This negligibly 
small effect of the canar.d on the wing loading 
due to wing pitching is In contrast with the 
significant effect of the canard on the steady 
wing loading (Fig. 13). Similarly, the 
oscillating wing has no effect on the fixed 
canard for separation distances of S = 0.50 and 
0.25 as shown In Fig. 17. This characteristic 
may be attributed to the relatively longer times 
for disturbances to propagate upstream from the 
oscillating wing to the motionless canard In 
comparison to downstream propagation from canard 
to wing as discussed previously for the 
wing-tall configuration. 


Concluding Remarks 

Steady and unsteady transonic flow 
calculations were presented for aerodynamlcally 
Interfering lifting surface configurations. The 
calculations were performed by extending the 
XTRAN3S unsteady transonic small -disturbance 
code to allow the treatment of an additional 
lifting surface. The code Is now capable of 
computing unsteady transonic flowflelds about 
relatively general wing-tall and canard-wing 
configurations. Grid generation procedures for 
swept and tapered Interfering lifting surface 
applications of XTRAN3S were described In 
detail. The planform finite-difference grid Is 
based on the single surface grid mapping 
procedures which were extended to Include an 
additional planform. Numerical Instabilities at 
the root sections which were encountered In 
Initial applications of XTRAN3S were eliminated 
using a new symmetric treatment of the grid 
across the reflection plane. Discontinuities In 
grid metrics were removed using an Inboard grid 
mapping procedure with parabolic smoothing. 

Unsteady subsonic calculations were 
performed for a simple wing-tall configuration 



Fig. 16 Unsteady sectional lift coefficients 
on canard-wing configuration due to 
canard pitching at M * 0.9, oq = 2°, 
and k * 0.3. 



Fig. 17 Unsteady sectional lift coefficients 
on canard-wing configuration due to 
wing pitching at M » 0.9, oq * 2°, 
and k * 0.3. 


10 



consisting of two identical swept planforms to 
assess the modifications to the XTRAN3S code. 
Unsteady chordwise pressure distributions and 
spanwise lift distributions for wing pitching 
showed good agreement with published linear 
theory results. Steady transonic calculations 
were performed for swept and tapered wing-tail 
and canard-wing configurations at several mean 
angles of attack. The effects of aerodynamic 
interference were demonstrated on transonic 
steady pressure distributions and spanwise lift 
distributions. For the wing-tail configuration, 
the wing produced a downwash on the tail which 
significantly decreased the tail spanwise lift 
distribution. Conversely, the tail induced an 
upwash on the wing which increased the wing 
spanwise lift. For the canard-wing 

configuration, the canard produced a downwash on 
the wing which decreased the Inboard wing 
spanwise lift at the higher angles of attack 
considered. The outboard wing lift was less 
affected by the presence of the canard, as 
expected. Unsteady transonic calculations were 
performed for wing, tall, or canard pitching 
motions. The effects of aerodynamic 

interference were demonstrated on transonic 
unsteady spanwise lift distributions. In 
general, the unsteady lift on one surface due to 
its own motion, termed self-induced lift, showed 
only small differences between Isolated and 
interfering lifting surface cases. The induced 
lift on the aft surface due. to motion of the 
forward surface, though, was dependent upon 
separation distance between the two lifting 
surfaces and was of similar magnitude to the 
self-induced lift. The effect of the 
oscillating aft surface on the forward surface 
was negligible. 

The research was conducted as a first-step 
toward developing the capability to treat a 
complete flight vehicle. Transonic aeroelastic 
calculations are presently being performed for a 
closely-coupled canard and forward-swept wing 
configuration to investigate the effects of 
aerodynamic Interference on stability and 
flutter. Work on developing a wing-body 
capability also is currently in progress. 
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